Directed polymers in random media under confining force 
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The scaling behavior of a directed polymer in a two-dimensional (2D) random potential under 
confining force is investigated. The energy of a polymer with configuration {y(x)} is given by 
H({y(x)}) = Y2x=i v{ x jV( x )) + eW a , where r)(x,y) is an uncorrelated random potential and W is 
the width of the polymer. Using an energy argument, it is conjectured that the radius of gyration 
R g (N) and the energy fluctuation AE(N) of the polymer of length N in the ground state increase 
as R g (N) ~ N v and AE(N) ~ N u respectively with v = 1/(1 + a) and u = (1 + 2a)/(4 + 4a) for 
a > 1/2. A novel algorithm of finding the exact ground state, with the effective time complexity of 
0(N 3 ), is introduced and used to confirm the conjecture numerically. 
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I. INTRODUCTION 



Scaling behaviors of a directed polymers in a random 
media (DPRM) have been studied extensively due to 
the mathematical interest of optimal path problems in 
random environment in addition to their application to 
physical systems such as stretched polymer in a gel or 
flux lines in a disordered superconductor [lj . A directed 
polymer is stretched in the longitudinal direction but can 
fluctuate in the transversal direction. We can map a con- 
figuration of a directed polymer in a ti-dimensional space 
to a path of random walkers in the d t — d — 1 dimensional 
space (in the transversal direction) when the longitudi- 
nal direction is considered as the time. Therefore, in suf- 
ficiently low temperature, a DPRM corresponds to the 
minimum energy path of a random walker in a random 
potential rj{t,r). Its scaling properties are characterized 
by two exponents, transverse length scale exponent v a 
and energy fluctuation exponent oj ||. For the d t — 1 
dimensional walkers, the exact values of these exponents, 
v = 2/3 and to = 1/3 can be obtained analytically by 
the combinatorial method |3J. The exact exponent val- 
ues can also be obtained by mapping the free energy of a 
random walker to Kardar-Parisi-Zhangequation (KPZ) 
via Cole-Hope transformation 0, 0, ■ This map- 
ping shows that the exponents for the DPRM, v Q and ui 
should be equal to the inverse of the dynamic exponent, 
l/z KPZ = 2/3 and the growth exponent, (3 KPZ = 1/3 
of KPZ equation respectively at any finite temperature. 
In other words, any disorder drives a directed polymer 
in the d — 2 dimensional media into a strong disorder, 
pinned phase at T = 0. 

In this paper, we study the scaling properties of a 
DPRM under confining force (DPRMCF). The Hamil- 
tonian for a DPRMCF has two terms, the usual random 
potential term Erm of the DPRM and the confining en- 
ergy term Eq which prefers the straight polymer in global 
length scale. It may describe a DPRM confined by an in- 
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FIG. 1: (Color online) Directed polymers in 2D random media 
under confining force. The magnitude of the random potential 
at each lattice site is denoted by the area of the circles at the 
lattice site. A directed polymer is confined by the two rods. 
The confining force may come from the elasticity which are 
symbolized by the spring (a) or from the repulsion of the 
polymers (red dashed lines) at the outside of the rods (b). 



flatable but non-flexible tube or by two rods with springs 
as illustrated in Fig. ^a) . We consider the confining en- 
ergy which depends only on the "global width" W of 
the polymer. Such confining energy term may mimic the 
elastic energy of the (inflatable) tube or the spring which 
prefers the smaller width. The confining force on the cen- 
ter polymer may arise from the repulsion from the other 
polymers (red dashed line) outside of the rods if we con- 
sider an array of polymers on a plate as illustrated in 
Fig. ^b). We may find other physical systems that our 
model might be applied, such as the motion of a single 
step on a terrace with quenched random impurities in a 
vicinal surface |8( but the principal motivation for the 
model is rather theoretical interest. The effects on the 
scaling properties of random walkers or directed poly- 
mers from the energy terms associated with the global 
configuration quantity, such as the global width W , are 
theoretically intriguing [{J, [lfj, [llj, [12J • For example, the 
confining energy term E c ~ W makes the normal random 
walkers (without Erm) visit the same y value sites even 
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times stochastically , and changes the "roughness" ex- 
ponent to 1/3 from the conventional universal value of 
1/2 d E El. The path of this stochastic even visit- 
ing random walker can be mapped a ID interface profile 
which is called a "self-flattening" surface |ToL Hll | . The 
height-height correlation function of the self-flattening 
surface shows an anomalous scaling behavior in the sense 
that the local wandering exponent is different from the 
global roughness exponent [ll| . This anomalous behavior 
indicates the existence of window length scale 1{N) ~ N s 
with < S < 1 above which the global energy term be- 
comes relevant. Recently, it has been conjectured that 



the window exponent S is given by 5 = 



D 



D+u 



from an 



energy-entropy argument where D is the dimension of the 
surface and v Q is the roughness exponent without the self- 
flattening mechanism |14j . This argument seems to be 
valid for general self-flattening mechanism, which corre- 
sponds to non-zero finite temperature dynamics with the 
confining energy proportional to the width W [TcHTHH^ . 

Here, we study the scaling properties of the zero tem- 
perature, ground state paths for a DPRMCF in 2D with 
a confining energy in the form of E c = eW a . If we as- 
sume the existence of window length l(N) ~ N s , a power 
counting argument predict S — 3/(2 + 2a) for a > 1/2 
and the confining energy of any positive e changes the 
roughness exponent v to 1/(1 + a) from v = 2/3 and 
the energy fluctuation exponent w to (1 + 2a) j (4 + 4a) 
from lo = 1/3. This conjecture is numerically tested for 
a series of different a values. 

This paper is organized as follows. In Sec. II, we de- 
fine the DPRMCF model. An analytic prediction for the 
roughness exponent v and the energy fluctuation expo- 
nent uj is presented in Sec. III. A novel algorithm for 
a numerical study of DPRMCF which finds the exact 
ground state in the effective time complexity of 0(N 3 ) is 
introduced in Sec. IV. In Sec. V, we present the numer- 
ical results from a series of Monte Carlo Simulations for 
DPRMCF model and calculate the exponents v and ui. 
We conclude the paper with a summary and remarks in 
Sec. VI. 



II. MODEL 

We consider a discrete model for a directed polymer in 
a random media under confining force (DPRMCF) whose 
Hamiltonian has two terms, site random potential and 
confining energy. Before introducing the details of our 
DPRMCF model, let us first explain a discrete model for 
DPRM without the confining force. We consider poly- 
mers on a 2D lattice whose horizontal axis is x and the 
vertical axis is y {i.e., the dt — 1 dimensional transversal 
space is the y-axis). At each lattice site, random po- 
tential rj(x, y) £ [0, 1] taken from a uniform distribution 
with 



(r)(x,y)T](x',y')) 



_L_ 

+ ^2 "x,x' Oy 



(1) 



is assigned. Since overhangs are forbidden for a directed 
polymer, its configuration is specified by a single valued 
integer function y(x) £ Z with |y(x+l) — y(x)\ — l. The 
energy of a DPRM is given by the sum of site potentials 
rj(x,y(x)) on which the polymer lies. At zero tempera- 
ture, a polymer on a given random potential locates at 
the path which minimizes the total site energies of the 
random media, 



N 



E RM ({y(x)}) = ^2n(x,y(x)), 



(2) 



where N is the length of the polymer. 

The configurational space of a DPRMCF is the same as 
that of a DPRM but the Hamiltonian H of a DPRMCF 
has the confining energy term Ec = e W a in addition to 
the Erm, 



N 



H({y(x)}) = J2r)(x,y(x))+eW a . 



(3) 



x = l 



Here, W is the absolute "width" of the polymer defined 
by 



W = y 



max 2/min 



i, 



(4) 



where y max and y m i n are the maximum and the minimum 
values of y(x) respectively and e and a are positive pa- 
rameters so that the confining energy increases when the 
width grows. 



III. WINDOW HYPOTHESIS AND POWER 
COUNTING ARGUMENT FROM ENERGY 
FLUCTUATION 

We study the scaling properties of the zero tempera- 
ture, ground state paths for the Hamiltonian of Eq. J2J . 
The statistical behaviors of a polymer in random media 
are usually characterized by the radius of gyration R g 
and the energy fluctuation AE defined by 



RgW = ] J((y-y) 2 ) 



1V 



\ in 



l N ( i N 

X — 1 \ X — l 



(5) 



and 
AE(N) 



E-(E) 



\ 



-y 



Hndv^xm-^H^ay^x)}) 



(6) 
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where f2 is the number of different realizations of random 
potentials and {y n (x)} and H v ({y n (x)}) are the ground 
state path and its energy for a given random potential 77. 
These two quantities asymptotically increase as 



v are given by 



R g (N) ~ N v 
AE{N) ~ A" 



(7) 
(8) 



for the ground state paths for a variety of different dis- 
tributions of the randomness 0, UiJ [l3 • 

In this section, we estimate v and ui for the DPRMCF 
by comparing the E R m and Ec in the ground state. For a 
directed polymer in a 2D random media, it has been well 
known that the energy of the minimum path of length N 
can be asymptotically written as 



E 



DPRM 



(N) =aN + bN u 



(9) 



with iV independent positive parameters, a, b and ui a . 
As N goes to infinity, the parameter a becomes sample 
independent and is about 0.25 for a random potential 
r](x, y) G [0, 1] taken from a uniform distribution. On the 
other hand, the second coefficient 6 is a sample dependent 
parameter with a finite variation. The energy fluctuation 
of the DPRM with length N is given by 



AEdprm{N) = \j\[E DPRM - (Edprm) j 

= ^{KbfN^, (10) 

and io is called the energy fluctuation exponent for the 
DPRM HEJ. 

We estimate v and ui for our DPRMCF model using 
Eq. with an assumption of the existence of "win- 
dow" |9l. ll l| . The confining energy term is assumed to be 
relevant only beyond a certain length scale of / ~ N s , so 
called "window" 0, El so that the polymer behaves like 
DPRM within the length scale of I. Then, the random 
potential energy E RM {1] N) of DPRMCF with length N 
can be written as 

E RM (l;N) = (al + br°)(N/l) 

= aN + bNl^"- 1 (11) 

and the total energy of DPRMCF is given by 

H(l;N) = E RM {N) + eW a 

= aN + bNl UJ °~ 1 +cl av °, (12) 

where a, b and c are independent of I and N since the 
absolute width of the polymer W is of order l v ° . From 
the minimization of H(l;N) with respect to I (i.e., from 
^ = 0), we have 

av 8-5 = l+u 5-25 (13) 

and the window exponent 5 and the roughness exponent 



S = 1/(1 + av Q - ujo) 

= 3/(2 + 2a), 

v = bv a 

= v j(\ + av Q - lu ) 



(14) 



= 1/(1 + a), (15) 

where we use v = 2/3 and lj q = 1/3 for ID DPRM @,H 
6] . The above two equation for S and v are only valid for 
a > 1/2 since 5 cannot be larger than 1. For < a < 1/2, 
we expect v — vq since the confining energy term is not 
relevant since the third term of Eq. l|12f) can be always 
neglected comparing to the second term as N goes to 
infinity. 

We can use a similar argument to estimate the energy 
fluctuation exponent. We first assume that the total en- 
ergy fluctuation AH should show the same N dependence 
as the random potential energy fluctuation AE R m for 
a > 1/2 from the power counting argument. Then, for 
the estimation of AE R m , we use the same "window" ar- 
gument used to extract the roughness exponent v. Since 
there are N/l ~ N x ~ s windows whose energies fluctuate 
with the "amplitude" l u ° ~ 7V Wo<5 , we have 



AE RM (N) = N u 



N- 



(16) 



and the energy fluctuation exponent uj is given by 

r 1 1 

UJ = LU 6 + - - - 

2 
l + 2a 



4(1 + a)' 



(17) 



IV. EXACT ENUMERATION ALGORITHM 

In this section, we consider a numerical method to 
find the ground state of the DPRMCF. For a given 
random energy r)(x,y) G [0,1] for the 2D lattice sites 
(x, y) £Z 2 , we find the minimum energy path of Eq. @ 
using an exact enumeration method under the constraint 
of \y(x+\) — y{x)\ — 1 with the anchored (y(Q) = 0) 
boundary condition. Since the number of all possible 
paths increases as 2 N for a polymer of length N, it is 
(computationally) impossible for large N (say, N > 50) 
to find the ground state by comparing the energies of all 
paths. We present a novel, efficient algorithm of finding 
the exact ground energy path with the effective time com- 
plexity 0(N 3 ), based on the transfer matrix algorithm 
(TMA) for the DPRM problem. Figure l^a) and (b) 
illustrate the TMA for the DPRM problem. We represent 
the random potential in (a) by writing its value at each 
lattice site in the circle at the site (For clarity, we show 
an integer-valued random potential rj(x, y) G {1,2,..., 9} 
instead of real numbers from [0, 1].) and show the mini- 
mum energy path in (b) by a solid line. The number in 
the circle at the site (x, y) in (b) is the minimum poten- 
tial energy E^(x, y) of a polymer from the origin to the 
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for \y\ < N and the ground state path is given by 



(8) \ 13 

S (9) i- 

\~~® 


15 




25 














FIG. 2: Directed polymer in 2D random media, (a) A random 
potential is assigned at each lattice site, (b) Transfer matrix 
algorithm for a DPRM. The number in the circle at the site 
(a;, y) is the sum of site potentials of the minimum energy 
path from the origin to the site (x,y). The solid line is the 
minimum energy path to x = N. (c) For the DPRMCF, the 
minimum energy path to the site (a;, y) is not a simple sum 
of the site (x,y) and the minimum energy path to (x—l,y—l) 
or [x— l,y+l) (see text), (d) The minimum energy path for a 
modified potential fj y y (x,y) with y L = —2 and y u = 1. 



site (x, y). The minimum energies E^ n {x,y) for all the 
(x, y) with \y\ < x for x G {1, 2, • • • , N} can be obtained 
in 0(N 2 ) time using the TMA, 

Em{x,v) = v{x,y) 

+ mm{E? n (x-l,y-l),E° m (x-l,y+l)} (18) 

with P^(0,0) = 77(0, 0) and E? n (x,y) = 00 for \y\ > x, 
where min[^4, B] is the minimum value of A and B. The 
minimum energy path P^(x,y) from the origin to the 
site (x, y) can be obtained similarly; 



(x,y)+P^(x-l,y-l) 

if E° m {x-l,y-l) < E° m {x-l,y+l) 
(x,y)+P^(x-l,y+l) 

if E° m (x-l,y-l) > E° m (x-l,y+l) 



(19) 



with P° (0,0) = (0,0). (If E° m (x-l,y-l) = E° m {x-l, y+1), 
the minimum energy path is not unique and both paths 
become the minimum energy paths, but this is a very 
rare event for real values of 77.) The ground state energy 
Eg(N) is then given by the minimum of E® n (N,y) over 



P°(iV) 



p 



{N,y° m (N)) 



(21) 



where y^(N) is the site such that E^(N, y^N)) = 
E g (N). The ground state path for the random poten- 
tial of Fig. |2f a) is obtained this way and shown in (b) as 
the solid line. 

One might suppose that a similar transfer matrix algo- 
rithm, with the time complexity of 0(N 2 ), might be pos- 
sible for a DPRMCF. The first term of Eq. ||THJl, r](x,y) 
is the energy cost for a DPRM to proceed one step fur- 
ther from x—1 to x. A naive generalization of Eq. 118|) 
is replacing r)(x, y) by r](x, y) + e [(W + l) a — W a ] when 
the one-step movement increases the polymer width from 
W to W + 1. However, such algorithm does not lead 
to the global ground state of Eq. J3J since the mini- 
mum energies E m (x— 1, y— 1) and E m (x—l,y+l) are not 
enough to determine E m (x,y) as illustrated in Fig. |5fc) 
where E m {x,y) is the energy (including the confining 
term) of the ground state path to the point (x, y) from 
the origin. As an example, consider the minimum en- 
ergy path to the point (4,2) P m (4, 2) for a = 1 and 
e = 2 with rj{x,y) given in Fig. ffia). By comparing 
energies of all paths, one can see that P m (4, 2) is given 
by P m (4,2) = (0,0)-(l,l)-(2,2)-(3,l)-(4,2), denoted by the 
solid line in Fig. |5fc) . This path is not a simple addition 
of the end point (4,2) to P m (3, 1) unlike the case of the 
DPRM. The minimum energy path to (3,1) is given by 
P ro (3,l)=(0,0)-(l,l)-(2,0)-(3,l), denoted by the dashed 
line. It has lower energy than P e (3,l) = (0,0)-(l,l)- 
(2,2)-(3,l) although 77(2,0) > 77(2, 2) since the global 
width of P m (3, 1) is smaller than that of P e (3, 1). How- 
ever, P m (4, 2) is given by the addition of the end point 
(4,2) to P e (3, 1) not to P m (3,l). This is because the 
global widths of both paths to (4,2), the path through 
P e (3, 1) and the path through P m (3,l), are the same 
as 3. Therefore, for a DPRMCF, we need information 
on the minimum height y m i n and the maximum height 
2/max of the path as well as the potential energy of the 
path. In other words, the minimum energy E m (x,y) up 
to the site {x, y) is not enough but we need to know 
the minimum energies E m (x,y,y m i n ,y max ) of the paths 
to (x, y) for all different combinations of y mm and y ma x- 
The minimum energy E m (x, y, y^, y max ) can be calcu- 
lated from the minimum energies E m (x—1, y-1, y' min , 2/max) 
and E m (x-l,y+l,y' min ,y' max ) for a proper combina- 
tions of y' min and y' maK values and the ground state en- 
ergy of Eq. ©, E g (N) is obtained as the minimum of 
E m (N,y,y min ,y max ), over y, y m [ n , and 



E g (N) 



mm 

V,Viain,Vn 



E m (N,y,y mini y max ). (22) 



E°(N) 



min E® n 

y 



(N,y) 



(20) 



In principle, one can construct an algorithm to find 
ground state energy and its path based on Eq. Il22[| but it 
generally requires 0(N 4 ) memory and time. Although its 
"effective" time complexity can be reduced 0(N 3 ) (see 
below), the 0(N 4 ) memory requirement puts a strong 
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upper bound on the sizes of the systems to be investi- 
gated. 

In this paper, we use a novel algorithm which uses only 
0(N 2 ) memory but is more efficient than the algorithm 
using Eq. (|22[) . As before, we imagine the set Y Vmin y m!a . 
of paths whose minimum and maximum heights are y m in 
and y max but we first consider the path with the mini- 
mum potential energy, 



mm 



N 



,x = l 



(23) 



instead of the minimum total energy. Since the minimum 
total energy £ m (A, y min , y max ) for the given y min and 
2/max values, is simply given by 

E m (N, 2/min,?/max) = #m (N, J/min, y max ) 

+ £ (ymax -ymin + 1)°, (24) 

the ground state energy E g (N) can be obtained by 

Eg (A) = min -^m(A, ymin, ymax) 

Smin.Smn 



mill [ E°(N, ymin , ymax) 

J/min,ymax 



Vu 



1) 



(25) 



from the minimum potential energies. Note that we al- 
ways have y m ; n < and y m ax > since we set y(0) = as 
the anchored boundary. Therefore, in general, we need 
to calculate E^(N, y m in, ymax) for A 2 different combina- 
tion of y min € {0, -1, . . . , —A} and y max G {0, 1, . . . , N} 
to get the ground state energy E g (N). However, for a 
confining energy with positive e and a, we do not have 
to look for the paths with y max > Wo or y mm < —Wo 
where Wo is the width of the minimum potential path 
Pg, which minimizes Erm of Eq. @. Let y min and 
2/max be the minimum and the maximum heights of the 
minimum potential path P°. Then E%(N, y min , y max ) is 
the minimum of the first term in the square bracket of 
Eq. © and Wo = y max - y lnin + 1. The total energy 
of a path, whose y ma x > Wo or y m - m < —Wo, cannot 
be smaller than that of Pg since its confining energy 
is larger than eW§ in addition to the fact that its po- 
tential energy is lager than £^(A,y° in ,y max ). There- 
fore, for a given r)(x,y), the ground state energy E g (N) 
of Eq. i125[) can be obtained from the minimization of 
E m (N, y min , y m ax) over y min e {0, -1, . . . , -Wq + 1} and 
y max e{0,l,...,Wo-l}. 

Now, let us introduce a simple way to calculate 

7 ymin i Umax 

) for the (y m in, ymax) pairs needed for 
the minimization of Eq. (|25|l . We use the conventional 
TMA for the DPRM problem but with a series of modi- 
fied site potentials, 



\x,y) 



r){x,y) for y L <y <y u , ^ 
oo otherwise, 



for y L > —Wo and y u < Wo- Note that we can- 
not obtained the minimum energies Pj^(A, y m i n , y max ) 
for all pairs of (y m in,yma X ) with y min > -Wq and 
ym ax < Wo by simply applying TMA to the poten- 
tial Vy L ,y u ^ smce (y m in,yma X ) of the minimum poten- 
tial path P® y of f]y L tVv is not necessary equal to 
(y L ,y r/ ). As an example, Fig. EJd) shows Py L . Vu for 
{VltVu) = ( — 2, 1). The minimum height of the Py L>y , 
ymin = — 2 is equal to y L but the maximum height of 
it, y m a X = is not equal to y u — 1. Therefore, the 
minimum path energies for the two modified random po- 
tentials ij-2,0 and 77-2.1 are the same as E® t (N, — 2,0) 
and we cannot obtain E? n (N, -2, 1) by the TMA with 
% ,j/ • However, this means that E^N, — 2, 1) is 
larger than £^(^,—2,0) and we can safely exclude 
E m (N, — 2,1) from the candidates of the E g (N) since 
E m (N, -2,1) must be larger than E m (N, -2, 0). In 
general, we cannot obtain E m (N,y L ,y u ) by the TMA 
with % y if (y r „in,y max ) of P° is not equal to 
(Ul ! 2/17 ) but we can exclude E m (A, y L ,y v ) in the min- 
imization of Eq. (|25l) since this means E m (N,y L ,y u ) is 
always larger than E m (N, y min , y max )- If, the (y min , y max ) 
value is not equal to (y t ,y u ), E^ l (N 7 y L7 y u ) is larger 
than ^(A, y min ,y max ) in addition to {y v -y L + 1)" > 
(y max - ymin + 1) Q - In other words, E g (N) can be ob- 
tained by finding the minimum of E m (N, y m in, 2/max) over 
only the (y m in,ym ax ) pairs which can be obtained by the 
TMA with fjy LtVu for y L > -Wq and y„ < Wq. 

Let us summary our algorithm to find the ground state 
for a DPRMCF. For a given r)(x,y), we first find the 
minimum potential path P° without the confining en- 
ergy term by using the conventional TMA for a DPRM 
and calculate Wo- This can be done in 0(N 2 ). Then 
we find Py^ y ^ by applying the TMA to a series of the 
modified potentials Vy L .y u for y L > —Wo and y v < Wq, 
and measure their potential energies E^ n {N, y m i n , y max ) 
and widths W = y max - y min + 1. For a given y L and y u , 
Py L , yu can be obtained in time of order N{y u —y L ). Since 

Wn ~ A 2 / 3 , all P° „ can be obtained in 0(N A ) on av- 
erage. The ground state energy E g (N) is then given by 
the minimization of E m (N, y m in, ym ax ) through Eq. I|25|) 
and the ground state path is given by the correspond- 
ing minimum energy path. Note that the ground state 
energies and their paths for all positive e and a values 
can be obtained using Eq. (|23|l once we get P° for 
y L > -Wo and y < Wq. 



V. NUMERICAL RESULTS 

We first consider a harmonic confining energy of e W 2 , 
that is, the a — 2 case with the random potentials 
r](x,y) £ [0,1] taken from the uniform distribution. We 
simulate a random potentials 77 using the computer- 
generated pseudo random numbers 18] and find the 
ground state path of Eq. using the algorithm pre- 
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FIG. 3: (a) Radius of gyration R g for iV=16, 23, 32, 45, 64, 
92, 128, 256, 362, and 512 systems with a harmonic confining 
energy E c = eW 2 (a = 2). From the top, R g for e=0, 0.001, 
0.002, 0.004, 0.008, 0.016 and 0.04 cases are shown. Fitting 
lines are in the form of R g = aN v with ^=0.65, 0.48, 0.43, 
0.38, 0.36, 0.34 and 0.33 respectively from the top. (b) Ef- 
fective roughness exponents for e=0, 0.001, 0.002, 0.004, 
0.008, 0.016 and 0.04 values are plotted against 1/N 1 with 
7 = 0.5 from the top. The uppermost curve of e = goes to 
the known value of v = 2/3 but all the other curves go to 
much smaller value around 1/3 as N goes to infinity. 



sented in the previous section. Then, we measure the 
square width R 2 (the square bracket of the Eq. JSJ) and 
the first E v and the second E 2 moments of the energy 
for the ground state path for the given rj. The radius of 

gyration R g and AE are then obtained as R g = (R 2 ) 

and AE = y/(E 2 ) - (E) 2 where (A) mean the average 
over different realization of random potentials. We use 
eight million different realizations (CI = 8 x 10 6 ) of ran- 
dom potentials for N=W, 23, 32, 45, 64, 92 and 128 sys- 
tems and four, two and one million different realizations 
for N=256, 362 and 512 systems respectively to obtain 
the average values. The average over this large number 
of different random potentials makes the statistical error 
bars smaller than the sizes of the symbols most cases ex- 
cept the effective energy fluctuation exponents oe shown 
in Fig. Efb) and Fig. \^d) later. 

FigureEJa) shows the radius of gyration R g for the har- 
monic confining energy E c = eW 2 as a function of poly- 
mer length N for e=0, 0.001, 0.002, 0.004, 0.008, 0.016, 



and 0.04 cases. For each e value, R g lies on a straight line 
in a log— log scale plot indicating R g (N) ~ N v . The least 
X 2 fits oiRg(N) ~ N v give ^=0.66, 0.48, 0.43, 0.38, 0.36, 
0.34 and 0.33 for e=0, 0.001, 0.002, 0.004, 0.008, 0.016, 
and 0.04 cases respectively. The "roughness" exponents 
for e > 0.04 are almost identical to those of the e = 0.04 
case unless e is very large where the finite size effect is 
strong. 

To estimate v values for N —> oo, we calculate ./V de- 
pendent effective roughness exponent v e ff defined by the 
successive slopes in the log— log plot. We use neighbor- 
ing three points to get the local slope, i.e., v e ff(Nk) is 
obtained as the slope of the least \ 2 fit using the three 
data Rg(Nk-i), R g (N k ), and R g (Nk+i), where N k are 
the system sizes in an ascending order, 16 = N± < N2 < 
■ ■ ■ < Nio = 512. The roughness exponent v is estimated 
by extrapolating u e ff in the infinite size limit with a plot 
v e ff against . In Fig. 01b), we choose 7 = 0.5 which 
characterizes the trend of v e ff for large N well. From the 
figure, we see that our numerical measurement of v for 
e = is consistent with the known value v = 2/3 of the 
DPRM. However, it is clear that the roughness exponents 
of the DPRMCF (e > 0) is not equal to v Q . It is difficult 
to extract the definite v values from our numerical data 
but we speculate that v e ff goes to the conjectured value 
of 1/3 for all e > as the system size N goes to infinity. 

Figure^a) shows the energy fluctuation AE as a func- 
tion of polymer length N for the above seven e values. 
As in the case of R g (N), for each e value, AE lies on a 
straight line in a log— log scale plot indicating a power 
law increase. The least x 2 fits of AE(N) ~ N u give 
w=0.32, 0.36, 0.37, 0.39, 0.40, 0.41 and 0.41 for e=0, 
0.001, 0.002, 0.004, 0.008, 0.016, and 0.04 cases respec- 
tively. We estimate the energy fluctuation exponent u> as 
the N — > 00 limit of the effective exponents oe as for the 
roughness exponent. The effective exponent oe is defined 
as the successive slopes in the log— log plot of Fig. 0Ja) 
and obtained by the least % 2 fit with the three neigh- 
boring points. In Fig. |I£b), the effective exponents are 
plotted against l/TV 7 with 7 = 0.5 as before. From the 
figure, we see that our numerical measurement of lu for 
e = is consistent with the known value lo = 1/3 of the 
DPRM. For e^0, all curves seem to go the conjectured 
value of u = ^ w 0.417. 

In addition to the radius of gyration, we also mea- 
sure the height (transverse position) correlation function 
G(x; N) to check the "window" argument. The corre- 
lation function G(x; N) is defined by the mean square 
height of the site a; of a polymer of length N, 

G(x;N) = (y(xf) 
{n(x,y)} 

where fl and y v (x) are the same as in the Eq. (j5j). 

FigurelSta) shows the G{x;N) for the DPRM (e = 0) 
for N = 32, 64, 128, 256 and 512 systems. All the data 
collapse to a single curve with G(x) ~ x 2 ^ 3 except the 
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FIG. 4: (a) Measured energy fluctuation AE for the systems 
of sizes AT=16, 23, 32, 45, 64, 92, 128, 256, 362, and 512 for 
a = 2. From the bottom, AE for e=0, 0.001, 0.002, 0.004, 
0.008, 0.016 and 0.04 cases are shown. Fitting lines are in 
the form of AE = aN" with o;=0.32, 0.36, 0.37, 0.39, 0.40, 
0.41 and 0.41 respectively from the bottom, (b) Effective 
exponents oe for the above 7 different values of e are plotted 
against 1/N 1 with 7 = 0.5. As N goes to infinity, the lowest 
curve of e = goes to the known value of ui = 1/3 but all the 
other curves go to much larger value around the predicted u> 
of 5/12. 



points near or at the boundaries at x = N. Note that 
the ground state path y n (x) for a given random potential 
77 cannot be determined locally even in the case of e = 
0. In other words, the minimum energy path up to x 
is different from the sub-path up to x of the minimum 
energy path of length N > x. Due to the global nature 
of the ground state path determination, G(x; N) deviate 
form the infinite size behavior for a finite portion of N 
even for the free boundary condition 0|. Yet, the height 
correlation function follows the scaling relation, 



G(x;N) = N 2 »° 9o (x/N) 



(28) 



where g (u) increases as g (u) ~ u 2v ° unless u is very 
close to 1 where the boundary effect exists. When we 
rescale the correlation function G(x; N) by N 2u ° with 
v Q = 2/3 and x by TV, all the data collapse to a single 
curve as shown in Fig. pf b) . 

However, the correlation function G(x; N) for e > 
shows qualitatively different behaviors from the DPRM 
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FIG. 5: (a) G(x) for the systems of sizes iV=32, 64, 128, 
256, and 512 with no confining force e = 0. (b) A rescaled 
height correlation G(x)/N 2v ° is plotted against rescaled dis- 
tance x/N with vq = 2/3 in log scale. 



case. When there is an energy term associated with the 
global width, there seems to be another length scale over 
which the correlation function saturated. Figure Eta) 
shows the G{x; N) for e = 0.04 for N = 32, 64, 128, 256 
and 512 systems. As x increases, the correlation func- 
tions increase algebraically only for x < l(N) and then 
remains as constant values. We rescale the correlation 
functions G(x;N) by N 2v and x by N s with the con- 
jectured values of Eq. (|15|) . v = 1/3 and 5=1/2 and 
plot G(x;N)/N 2v against x/N s in Fig. Et*>). All data 
collapse to a single curve, implying a new scaling law 

GO; TV) = N 2u g{x/N s ). (29) 

The scaling function g(u) increases as g(u) ~ u 2v ° for 
u < 1 and then becomes a constant for u > 1. In other 
words, G(x; N) increase as x 2v ° for x < N s and reaches a 
constant value for x > N s . Note that the scaling function 
grows algebraically with exponent 2^o not with 2u sup- 
porting the assumption that our polymer behaviors like 
a DPRM up to the window size N s and then feels the 
global constraints of the confining force over the window 
size. 

We perform a series of simulations for other values of 
a and obtain the numerically results consistent with the 
conjecture of Eq. (|15fl . Figure [7| shows the radius of gy- 
ration R g and the energy fluctuation AE and their ef- 
fective exponents v e ff and ce for the confining energy of 
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FIG. 6: (a) The correlation function G(x) for the systems 
of sizes iV=32, 64, 128, 256, and 512 with a confining energy 
E c = 0.04W 2 . (b) The rescaled height correlation G(x)/N 2v 
is plotted against rescaled distance x/N s in a log scale with 
v = 1/3 and 5 = 1/2. 



W c = eW (a — 1). The radius of gyration R g lies on a 
straight line in a log— log scale plot as before and the least 
X 2 fits of R g (N) ~ N u give ^=0.65, 0.57, 0.53, 0.49, 0.46 
and 0.45 for e=0, 0.05, 0.1, 0.2, 0.5, and 1.0 respectively. 
Note that the v values from the least % 2 fits (using R g 
for N < 512) are smaller than the conjectured value of 
1/2 for some large epsilon values while those for small 
cpsilon values are larger than the conjectured values un- 
like the case of a = 2. Yet, as the system size goes to 
infinity, all the effective roughness exponents v e ff seem to 
go to the conjectured value of 1/2 for e > while i/ e ff for 
e = goes to the known value of 2/3 (see Fig. EJb)). The 
energy fluctuation AE is also measured and the effec- 
tive energy fluctuation exponents ce are calculated from 
them. As shown in Fig. [7f c) , the energy fluctuation AE 
increases as AE ~ N u with w=0.32, 0.33, 0.34, 0.35, 0.35 
and 0.35 for e=0, 0.05, 0.1, 0.2, 0.5, and 1.0 respectively. 
From these data, it seems to be difficult to distinguish 
the energy fluctuation exponent of DPRMCF from that 
of DPRM. The analysis of the effective exponents v e ff(N) 
provides somewhat better diagnosis. In Fig.[7Jd) we plot 
v e ff against 1/N - 5 as before. The error bars of ce for 
a = 1 are relatively large as shown in the figure although 
we obtain the effective exponent with more than million 
different random potentials. The statistical error of one 
millionth in measuring the second moment average (E 2 ) 




FIG 



for the systems of sizes N- 



32, 45, 64, 92, 128, 256, 362, and 512 for a = 1. From the 
top, R g for e=0, 0.05, 0.1, 0.2, 0.5, and 1.0 cases are shown. 
Fitting lines are in the form of R g — aN v with !/=0.65, 0.57, 
0.53, 0.49, 0.46 and 0.45 from the top. (b) Effective roughness 
exponents for the above 7 different values of e are plotted 
against 1/iV 7 with y = 0.5. As the system size goes to in- 
finity, the effective exponents seems to go to the conjectured 
value 0.5 for all e / 0. (c) Energy fluctuations AE are plotted 
against N. Fitting lines are in the form of AE = aN^ with 
w=0.32, 0.33, 0.34, 0.35, 0.35 and 0.35 in an ascending order 
or epsilon values from the bottom, (d) Effective energy fluc- 
tuation exponents oe against 1/N 1 with 7 = 0.5. The lowest 
curve is for the e = case. 



gives rise to the error bars of the symbol size in Fig. Eld). 
Due to these limitations, we cannot extract the definite 
value of u) for e > from our simulations but the nu- 
merical results does not seem to exclude the conjectured 
value of u) = 3/8 = 0.375. Our simulations on other a 
values such as a — 1.5 and a = 3 also give the consistent 
results with the conjecture of Eqs. I|15f) and (|17fl . 



VI. CONCLUDING REMARKS 

We consider the scaling behavior of a directed polymer 
in a 2D random media with confining energy E c = eW a 
and find that the roughness exponent v and the energy 
fluctuation exponent u> are given by v = l/(a + 1) and 
u> = (2a + l)/4(a + 1) respectively. These results can 
be understood by assuming that a polymer of length N 
behaves like a DPRM up to the window size 1{N) ~ N s 
and then feels the confining energy over the window size. 

We have only considered the scaling behavior of the 
zero temperature ground state polymers. We know the fi- 
nite temperature behaviors of the polymers only for some 
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limiting cases. For DPRM where confining energy term is 
absent, zero temperature pinned phase is the fixed point 
so that polymers at any finite temperature shows the 
same exponents with those at the zero temperature. On 
the other hands, if there is only confining energy (without 
random potential) there are three phases as temperature 
changes. At zero temperature, the polymer becomes a 
straight line with width 1 and therefore v — while it 
becomes a random walk with v = 1/2 at infinite tempera- 
ture. At nonzero finite temperature, especially for a = 1 
where the confining energy is given by E c — eW, the 
polymer configurations are identical to the self flattening 
surface whose roughness exponent v is 1/3. When 
there are both random potential and confining energy 
terms, we do not know the scaling behaviors of the poly- 
mers at finite temperature where our algorithm cannot 



be applied. Further investigations are needed to explore 
full phase diagram of DPRMCF over general temperature 
and dimensions. 
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